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Biological organisms swimming at low Reynolds number are often influenced by the presence 
of rigid boundaries and soft interfaces. In this paper we present an analysis of locomotion near 
f^*) , a free surface with surface tension. Using a simplified two-dimensional singularity model, and 

£\j ■ combining a complex variable approach with conformal mapping techniques, we demonstrate that 

the deformation of a free surface can be harnessed to produce steady locomotion parallel to the 
interface. The crucial physical ingredient lies in the nonlinear hydrodynamic coupling between the 
disturbance flow created by the swimmer and the free boundary problem at the fluid surface. 
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(N| ■ I. INTRODUCTION 

Low-Reynolds number swimming near solid boundaries or interfaces can exhibit interesting and unexpected features. 



In particular, the presence of long-range interactions typical of flows at low Reynolds numbers implies that, in general, 
boundary effects can not be ignored For instance, E. coli cells arc observed to change their swimming trajectories 

from straight to circular when they are moving parallel to a solid surface f^-IB], a behaviour modification which may 
have important implications in the formation of biofilms The motion of microorganisms near soft interfaces, 



i 

q 

^ ■ such as spermatozoa motility through the mucus- filled female reproductive track [8j, is even more intriguing as the 
* J5 ■ nonlinear coupling between the motion of the swimmer and the changing shape of the interfaces adds an extra level 
of complexity to the problem. 

The study of low Reynolds number swimmers near a no-slip wall has received considerable attention in the past, 
1—1 . and we refer to the reviews by []J and Q for a discussion of the relevant literature. Most theoretical work has 
, focused on quantifying the change in swimming speed and energetics near solid boundaries More recent work 



has addressed the dynamics of confined swimmers, and tackled the subtle interplay between the time evolution of a 
swimmer's orientation and its position. For example, a well-known feature of swimming near a solid boundary is that 
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t^J- . organisms moving at low Reynolds number tend to be attracted to solid surfaces jl3l - ll9| . This phenomenon can be 
rationalized by a fundamentally hydrodynamical mechanism in which the interaction with the rigid boundary causes 
a swimmer to reorient itself in such a way that it is eventually attracted to its hydrodynamic image system in the 

§ '■ wall 0- 

Other studies have revealed additional dynamical features of a swimmer's behaviour near a wall, f^of have conducted 
numerical experiments to understand the wall-bounded dynamics of model swimmers from a control and dynamical 
systems perspective. In addition to the existence of a steady state in which the swimmers travel in a steady rectilinear 
motion parallel to the wall, the authors found that the generic motion of a swimmer can be described by nonlinear 



* 

. . i periodic orbits along the wall with complicated spatio-temporal structure. These observations have since been cor- 
roborated by laboratory experiments involving small robotic swimmers in a tank of viscous fluid (2l| . Motivated by 



these studies, 22( have recently proposed a simple two dimensional model of a swimmer near a wall. Using a complex 
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variable formulation of the problem, they obtained results in agreement with those of [2fJ and 

In the current paper, we address the coupling between a low-Reynolds swimmer and a surface which can deform, 
and focus on the case of a free interface with surface tension. Previous work considered how the unsteady deformation 



of soft surfaces generated by time-reversible flows could provide new modes of locomotion and pumping [23[. Here 
we ask the following question: Can a low-Reynolds number body exploit the deformation of a free surface to swim 
steadily? To emphasize the role played by the surface deformation, we consider swimmers which cannot swim in the 
absence of a free surface, and determine whether they are able, through the disturbance flow field they are creating 
and the subsequent surface deformation, to acquire locomotive abilities. Specifically, using the modeling approach of 
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22 , we focus on identifying a steady mode of locomotion in which a swimmer translates at a constant speed parallel 



to the undisturbed free surface. Since, in the neighborhood of a flat no-slip wall, the motion of a swimmer generically 
follows a time-dependent periodic orbit, it is not clear a priori that a steadily translating swimmer motion beneath 
a free capillary surface is possible. What we show below is that such a mode of locomotion is indeed possible and, 
within our model can be described in a mathematically explicit way. 

Our study was originally motivated by the discovery of a peculiar mode of locomotion employed by water snails 
that crawl underneath the free surface. Separated from the interface by a thin layer of mucus, these organisms deform 
their foot to create a lubrication flow inside the mucus layer. This flow results in deformations of the free surface, 



which in turn rectifies the flow, allowing the water snails to move [241 ] . The analysis in this paper contained two 
significant constraints, namely the gap between the swimmer and the air-water interface was assumed to be thin, and 
the deformation of the free surface was assumed to asymptotically small. The work in the current paper removes these 
constraints and considers a more general mechanism for a swimmer translating steadily beneath the free surface. 

Low Reynolds number swimmers exert no net force and no net torque on the flow, and it is precisely these constraints 
that dictate the subsequent speed of the swimmer and its angular velocity. Here we introduce a mathematical 
representation of the swimmer as a two-dimensional torque-free point stresslet which, by definition, is force and 



torque-free. This type of singularity model has been widely used in modeling suspensions of force-free particles [25| 



and swimming microorganisms 13, l26l . 1271 ] . The approach is equivalent to considering the swimmer on distances much 



larger than its intrinsic size, so that its precise geometric structure and the fine details of its swimming protocol are 
encapsulated in the effective far-field multipole structure. The two-dimensional assumption, although idealized and 
not directly relevant to biological swimmers, allows us to explicitly solve for the nonlinear free boundary problem, 
thereby shedding light on this new mode of locomotion. 

Our mathematical approach is inspired by the work of 28j] who considered surface deformations generated by 
two counter-rotating cylinders beneath a free surface at low Reynolds numbers. In a similarly idealized model, the 
flow generated by the two counter-rotating cylinders is modeled by a single potential dipolc located on the axis 
of symmetry of the deformed free surface. This flow results in symmetric deformations of the interface which are 
calculated, for a given dipole strength and fluid properties, by means of conformal mapping methods. Notably, the 
conformal map approach allowed them to produce exact solutions even for large nonlinear deformations of the free 
surface. Furthermore, their analytic results exhibit remarkable agreement with the experimental data of free surface 
deformations for different rotation rates of the cylinders. In the spirit of 28( we therefore use a conformal map to 
solve for the shapes of the free surface governed by the interaction between surface tension and the flow field generated 
by the swimmer, represented as a combination of singularities. The difference here is that we must allow for non 
symmetric deformations of the interface and adapt the analysis to admit a stresslet singularity in the fluid (rather 
than a potential dipole). 

This paper is organized as follows. In §[TT1 the two dimensional Stokes equations and relevant boundary conditions 
are introduced in complex variables. The singularity model approach is then explained in § 11111 Section IIVI uses 
a method of images to demonstrate that steady motion of a swimmer beneath a flat undeformed interface is not 
possible. Section [V] then introduces a conformal mapping approach that enables us to explore solutions in which 
the free surface admits essentially arbitrary deformations. Section IVII gives a characterization of the class of steadily 
translating solutions, and is followed by the conclusions in ij lVIII 



II. COMPLEX VARIABLE FORMULATION OF STOKES FLOW 



Let the two-dimensional quiescent fluid occupy the area beneath a deformable fluid-air interface, D. The fluid is 
assumed to be incompressible and, in the Stokes regime, the streamfunction ijj is known to satisfy the biharmonic 
equation 

V 4 i>(x,y) = 0. (1) 
Introducing the complex- valued coordinate z = x + vij, it is possible to write the general solution of the biharmonic 



3 



equation in the form 

4> = Im{!f(z)+g(z)}. (2) 

Here / = /(£) and g = g(z) are two functions which must be analytic functions of z inside the fluid region except 
at isolated points where singularities arc deliberately introduced in order to model particular flow conditions. These 
functions are sometimes referred to as Goursat functions. 

It is possible [2^ to express all the usual physical variables in terms of these two functions. Indeed, it can be shown 
that 

P 



iu> = if'(z), 



u + w = -f{z) + zfi(z)+g>{z), 



(3) 



e n +ie 1 2 = zf"{z)+g"{z). 



Here, p is the fluid pressure, ui is the vorticity, (u, v) is the fluid velocity and is the fluid rate-of-strain tensor. The 
dynamic fluid viscosity is p. Primes denote differentiation with respect to z, and overbars denote complex conjugates. 

The stress boundary condition on the free surface requires that the normal fluid stress is balanced by the surface 
tension and that the tangential stress vanishes. This can be written as 

- prii + 1'/"' <.;». = (4) 

where a is the surface tension, k is the surface curvature, and rii is the outward unit vector normal to the interface. In 
addition, the kinematic condition on the interface requires that the normal velocity of the interface equals the normal 
fluid velocity. 

The governing equations, Goursat functions, and corresponding boundary conditions are non-dimensionalized as 
follows: 

z = hz, z c i = hz c i, u + iv = U(u + iv), -0 = Uhip, 

- . aU ( 5 ) 

f = Uf, g= Uhg, p = -j-p, 

h 

where z^ is the dimensional location of the swimmer, h is the magnitude of Zd or the vertical distance of the swimmer 
from the interface, and U is a characteristic speed of translation. The capillary number Ca, which reflects the 
dimcnsionlcss ratio of viscous to capillary effects, is defined as 

Ca=^. (6) 
a 

It can be shown that the complex form of the stress condition (j?]) on the air-fluid interface is equivalent to the relation 

dH i d 2 z 

where ds is a differential element of arc length along the free surface and 

H^f(z) + zj\z)+g'{z). (8) 
Hence, the stress condition can be integrated once with respect to s to give 

f{t) + z ?W+ m = -^*L, (9) 



where, without loss of generality, the constant of integration has been set equal to zero. 
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FIG. 1: Illustration of the singularity model: A finite-body swimmer beneath a free surface is modelled as a point stresslet 
singularity with superposed potential dipole and quadrupole. 



III. A SINGULARITY MODEL OF THE SWIMMER 

From the first equation in §3§ it is clear that singularities of the Goursat function f(z) will be related to local sin- 
gularities in the pressure field and hence, to localized force singularities. A logarithmic singularity of f{z) corresponds 



to what is often referred to as a stokeslet, or point force singularity 30( which, as previously discussed, is not allowed 
owing to the condition that the swimmer exerts no net force on the flow. 

The next order singularity (the derivative of the logarithm) is a simple pole. If near Zd, f(z) has a simple pole 
singularity, 

s* 

f(z) = h analytic, (10) 

z - z d 

then, in order to ensure that the velocity field scales like l/\z — z d \ (rather than l/\z — z d \ 2 ) we must also have 

g'{z) = — — — ^— - + analytic. (11) 

[z z d ) 

Thus, if f(z) and g'{z) locally have the behavior (fT0|) and (JTTJ) respectively near Zd then there is a stresslet of strength 
s* at Zd- In general, if g(z) has a simple pole near some point Zd of the form 

g(z) = — h analytic, (12) 

z - z d 

then we say there is a dipole of strength d at Zd- A stresslet singularity of strength s* at Zd therefore corresponds to 
a simple pole of f(z) at Zd with residue s* together with a simple pole of g{z) with residue — s*z~d at the same point. 

There is a physical way to understand this singularity of f(z). Swimmers at low Reynolds numbers propel themselves 
by exerting a local force on the flow (for example, the waving flagcllum of a spermatozoa) which is then counterbalanced 
by a net drag on its body (the head of the spermatozoa plus its flagcllum) leaving the total net force on it equal to zero 
[e.g. see review inQ. If this scenario is modeled as two logarithmic singularities of f(z) (two point forces) drawing 
infinitesimally close together with equal and opposite strengths tending to infinity at a rate inversely proportional to 
their separation, the limit is precisely a simple pole singularity of f(z) of the form (|10p . 

In a general singularity description of a swimmer, g'(z) is also singular at Zd- We already know it must have a 
second order pole (fTTj) (this is associated with the stresslet) but it can have additional singularities. The singularities 
of g(z) are potential multipolcs because, as is clear from (j3|), only the singularities of f(z) contribute to the vorticity of 
the flow. Different swimmers generate different effective singularities according to their particular swimming protocol. 
Any choice of singularities that we assume g'(z) to have at z d is therefore a manifestation of our choice of swimmer 
type. It is not clear, a priori, how to pick either the type of these singularities of g'{z) or their magnitudes. 

In this paper wc adopt the same singularity model of a swimmer used by [2^] in their studies of a low-Reynolds 
number swimmer near a no-slip wall. They motivated their choice of singularities by considering a concrete model of 



5 



a finite-area circular "trcadmilling" swimmer of radius e. It was supposed that, on its surface, the swimmer generates 
a purely tangential surface velocity given by 

U(</),t) = 2Vsm(2(<f>-6(t))) (13) 

where V is a constant (setting the time-scale of the swimmer's motion), <f> is the angular variable and 9(t) is a 
distinguished angle taken to be the direction in which the head of the swimmer is pointed. By solving a boundary 
value problem for the flow associated with this swimmer in an unbounded Stokes flow, it is possible to show that such 
a swimmer has an effective singularity description consisting of a stresslet of strength fi(t) = eV exp(2i#(f)) with a 
superposed potential quadrupole of strength 2^{t)e 2 . 

This model is a particular case of a general class of simplified swimmers first considered in a theoretical study due 
to who looked at the effect of imposing velocity profiles of general form on the surface of a circular swimmer. Such 
an "envelope model" captures the macroscopic effect of the motion of many small-scale beating cilia on the swimmer 
surface. Similarly, cilia-aided crawling of organisms beneath a free surface has been observed in nature, in particular 
for some families of snails. [32J concluded that the locomotion of Alectrion trivittata, which crawls upside down on the 
surface, relies solely on the ciliary action. He conducted a similar study on Polinices duplicate/, and Polinices heros, 
both of which were observed to use both cilia and muscle contraction for locomotion on hard surfaces 33j. Only 



ciliary motion was employed by the young Polinices heros when crawling inverted beneath the surface. 

Prompted by the success of the previously described singularity model, we extend the study to point swimmers 
(beneath a free surface) within the same general class: that is, a point stresslet superposed with a potential dipole 
and quadrupole. The dipole has been included because it is a lower order singularity than the quadrupole and there 
is no reason a priori to suppose it is absent (moreover steady rectilinear motion is expected to involve a dipole in its 
singularity description). This means that we will seek f(z) and g(z) with the functional forms 

f(z) = — - + fo + fi(z + i) + ..., g(z)= q + , +g + gi (z + i) + ... (14) 

z + i [z + iy (z + i) 

where the singularity is at z = — i and /o, /i, go and g\ are constants. We will refer to s* as the stresslet strength, q* as 
the quadrupole strength and d* as the dipole strength (note that part of the coefficient of 1/(2 + i) in g(z) is naturally 
associated with the stresslet singularity as seen in (|T0|) and ([XT])). We will not make any a priori assumptions on the 
relative magnitudes of s*,d* and q* since, for a steady solution, we expect these to be determined by the conditions 
fo r eq uilibrium. 



22j show that the evolution equations for the swimmer position Zd(t) and its orientation 9(t) are given by the 



dynamical system 



d -^ = - h + Zd J l+Tl , « = -2Im[/ l] . (15) 



The first equation states that the swimmer moves with the finite part of the fluid velocity at the swimmer position, 
while the second equation states that its angular velocity equals half the regular part of the vorticity at the swimmer 
position. In the present paper we adopt these same evolution equations but focus on finding equilibrium solutions in 
which the swimmer translates steadily in the direction of the undeformcd interface (i.e., parallel to the a;-axis). In a 
co-travelling frame we therefore need to find solutions satisfying the conditions 

= -fa + Zih + W, = Im[/i] . (16) 

The first equation ensures that the swimmer is stationary in the co-moving frame. The second equation ensures that 
the local vorticity at the swimmer position vanishes so that its orientation remains fixed in time. Note that for a 
no-slip wall the time evolution of the swimmer's orientation proves to be a crucial ingredient in understanding the 
dynamics of swimmers near those surfaces |l9| , and the same is expected to be true near a free surface. 
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FIG. 2: The method of images near a flat free surface: to satisfy the boundary condition on the surface an image singularity 
must be introduced at the reflection of the point swimmer position. 



IV. STEADY SWIMMING BENEATH A FLAT INTERFACE: METHOD OF IMAGES 



It is instructive to first examine whether it is possible to find a solution for a point swimmer (of the type just 
described) translating steadily beneath a flat undeformed interface. To do so we must seek /(z) and g(z) satisfying 
the boundary conditions on the free surface, with singularities of the form (|14[) . and satisfying conditions (|16[) . Owing 
to the simplicity of the geometry we search for such a solution using the method of images, as originally introduced by 
Blake and co-workers for three-dimensional flows near no- slip walls 0,0. We now demonstrate that such a solution 
does not exist. In subsequent sections we succeed in identifying steadily translating solutions in which the interface 
deforms. This highlights the crucial role played by the deformability of the interface in providing a mechanism for 
steady translation of the swimmer. 

Consider a swimmer moving at a prescribed horizontal speed U = — 1 beneath a flat interface in a Stokes flow with 
no background flow. It is convenient to work in a frame traveling with the swimmer. Then the swimmer is stationary 
at Zd = — i while the fluid in the far-field has unit speed in the x-dircction. 

We seek a solution of the form 



/(*) 



z + i 



- + /o 



(17) 



where, in addition to the stresslet of strength s* at zi = — i, we have placed an image stresslet of strength s* (to be 
determined) at z = i and fo is a constant. Figure [2] shows a schematic of both the swimmer and its image. The form 
of g'(z) is now forced by the stress boundary condition on this interface. On z = zwe have, from ([9jl. 



g'(z) = -f(z) - zf (z) + 



2Ca' 



(18) 



Substitution of (fl~7f into (fl"8f and picking s* = s* so that g'(z) has no rotlet contribution (no simple poles) leads to 

is* is* 



9\z) = 



(z + i) 2 (z-i) 2 



,9o 



(19) 



which we rewrite as 



(z + i) 2 (z-i) 2 



2is* 



2is* 



(z + i) 2 (z-i) 2 



5o- 



(20) 



As can be seen by comparison of (|T7|) and (|2T)|) with (fTU|) and (|TT|) . the two second order poles in square brackets 
are associated with the two stresslets (the swimmer stresslet and its image) and the additional second order poles 
correspond to superposed potential dipole singularities (one at the swimmer position and another at the image point). 

The constants fo and go must be picked so that the far-field velocity condition is satisfied. Making use of (flT)) and 
(|T5)) in the expression for the velocity §3§ in the far field equation produces 



fo + .9o = 1. 



(21) 
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By balancing the constant terms in the stress condition, we can deduce 

/o = ~2 + 4(b go = \~lk- (22) 
Finally, the swimmer must be stationary in the cotranslating frame. Hence the finite part of 

-f( z ) + zf'(z)+g'(z) (23) 

at z = — i must vanish. This leads to 

°=-(^- /o+ (5r(5) + ^ (24) 

and the conclusion that 

s* = 2i. (25) 

Note that the strength of these singularities does not depend on the capillary number. 

The solutions for f(z) and g'(z) just derived give the instantaneous velocity field. We must also check that the free 
surface is a streamline. It can be verified (see the discussion in the next section) that this condition is equivalent to 

g(z)+zf(z)=0 (26) 

on the real axis. By making use of (fT7|) and ([T8|) we find 

is is* zs zs* 

g(z) + zf{z) = - — — - + zg + + + zf 

(z + i) (z-i) (z + i) (z-i) 

= 2Re[s*] = 

which confirms that the free surface is indeed a streamline. 

It still remains to satisfy the condition that Im[/i] = but it is easily verified that this condition is not satisfied. 
Furthermore, there are no further degrees of freedom in the solution that will allow us to enforce it. We therefore 
conclude that steady, non-rotating solutions for a steadily translating swimmer beneath an undeformed interface do 
not exist. 



V. STEADY SWIMMING BENEATH A DEFORMED INTERFACE: CONFORMAL MAPPING 

The failure to identify a solution for a steadily translating swimmer beneath an undeformed interface does not 
preclude the existence of such solutions when the interface deforms. We now seek such solutions. When the interface 
is not flat it is no longer obvious how to apply the method of images. Instead, we introduce a new solution method 
based on conformal mapping. This enables us to successfully resolve the nonlinear interaction between the point 
swimmer and the deformed free surface. 

Consider the unit disc |£| < 1 in a complex parametric £-plane. It is known, by the Riemann mapping theorem [3^), 
that there exists a conformal mapping z(Q that will map this unit disc to the fluid region beneath the free surface. 
The free surface itself will be the image of the unit circle \(\ = 1 under this mapping (see Figured]). Since the physical 
fluid domain is unbounded and the free surface extends to infinity, there must be a simple pole of z(£) on the unit 
^-circle. Using a rotational degree of freedom of the Riemann mapping theorem, we take this pole to be at £ = — i. 
The most general form of the mapping can then be represented as 

00 

^C) = -^T + £ a ^ fe (27) 

^ 1 k=Q 
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FIG. 3: Conformal mapping z(() from the unit ("-disc to the fluid region beneath the interface. In this mapping, £ = maps 
to the swimmer at z = Zd, and f = — i maps to the interface at infinity. The variable s denotes the arclength and increases as 
shown in the figure. 



where a and {a^} are a (generally infinite) set of complex coefficients. For a one-to-one function, it must also be true 
that dz/dC, ^ inside the unit disc. There are two remaining degrees of freedom in the mapping theorem which allow 
us to prescribe that the singularity Zd corresponding to the swimmer is the image of £ = 0, i.e, 



z d = z(0). 



(28) 



This means that the pre- image of Zd is the point £ — 0. 

We seek a solution where the singularity travels uniformly with speed U = — 1 in the i-direction and we again move 
to a co-travelling frame in which both the swimmer and the shape of the free surface are stationary. The kinematic 
condition on the interface for a steady solution in this frame is 



u n = 0. 



(29) 



If the normal vector n = (n x , n y ), the complex normal n x + \n y is —idz/ds, where s increases along the interface from 
positive infinity to negative (see Figure [3]). Using this fact, (|2l?)) can be expressed in the complex form as 



Re 



, . ,.dz 
(u + w li- 
as 



0. 



(30) 



This condition is equivalent to the free surface being a streamline and, in turn, corresponds to the following condition 
at the free surface: 



zf(z) + g(z) = 0. 



(31) 



The easiest way to show this is to compute the derivative (with respect to arclength s) of the quantity zf(z) + g(z) 
on the interface and then make use of © and (pttlf to show that ~zf(z) + g(z) equals zero. 
In the far-field it is required that 

(!) „(1) 



and 



g'(z) -> Soo + ^ + . . . 



(32) 



so that, as z — > oo, 



u + iv ^ -f^+g^+Odz^ 1 ). 



(33) 



To satisfy the boundary conditions at infinity that the interface moves in the cc-direction with unit speed we must 
have 



foo + g a 



i. 



(34) 



It can be shown that 



A + lk goo = l-JL-. 



2 4Ca 

To see how these relations arise, we make use of the fact that 

dz i£z£ 
ds \z{\ 



and show that, as C — > ~h 



so that, using the condition (J5), as £ — >• — i, 



Since, as ( — ► — i on \(\ = 1, 



i dz i a 



2 ds 2Ca a 



foo + 9o 



i a 
2Ca ]a\ 



- -f/—i\ a C , — , "i a 



ao 



_ , , ai 



1-iC C C + i 



t + la + a + 



it also follows, from ([3T]). that 



and hence 



This means, from that 

_ r _ /oO fl _ -, 

Joo _ — 1- 

a 

The only way for the velocity to be purely real is for a to be real so that 

foo foo 1 • 

It further follows, from (1381). that 



foo foo 2Ca' 

leading to (|3"5j) . 

Together with the stress condition ([9]), the kinematic condition (|29|) can be written 



Rc 



, dz 

2 ^dS 



l 

2Ca' 



We now introduce the composed functions 

F(Q = f(z(Q), G(0=9(z(0), 
which can be used in the kinematic condition together with (|36f to give 



Re 



' mo ' 



2Ca|z c | ; 



10 



where 2^(C) = dz/dC,. Since f(z) is required to have a simple pole at z^, F(C) must have a simple pole at £ = 0. 
Therefore 



F(Q = ^ + F + O(0 



for some constant Fa (to be determined). Therefore, consider 

1 



Re 



F(Q D C 
C*c C 2 C 



4Ca|z f | 



Re 



D C 

c 2 + c 



4Ce 



Re 



DC 2 + C( 



(48) 



(49) 



For suitable choices of constants C and D the function in the square brackets on the left hand side of (|49| is analytic 
everywhere inside the unit C-disc. Equation (|49p therefore gives the real part of an analytic function on the boundary 
of the disc. The Poisson integral formula can be used to give 



F(C) = (z c (C) 7(£,Ca) 
where b is some real constant and 

7(C,Ca) = 



D C 



DC ~C( + ib 



< C' + C i 
87riCa 7j c ,| =1 c C'-CWCOr 



(50) 



(51) 



Since z$ (C) has a second-order pole at ( — — i, it is necessary that the quantity inside the brackets in expression (|50l) . 
and its derivative, vanish at £ = — i so that near £ = — i 



7(C, Ca) + § + y - £C 2 - C( + ib = A(( + i) 2 + 

c c 

for some complex constant A which is related to /oo by 

iaA = foo. 

This means that 

7(-i, Ca) - D + D + iC + iC + ib = 0, 

and 

7 c (-i, Ca) + 2iD + C + 2Di - C = 0, 

while 

7 Cf (-i, Ca) + 6D- 2Ci - 2D 



A='- 
2 



(52) 



(53) 



(54) 



(55) 



(56) 



Since 7(— i, Ca) is purely imaginary, (|54[) is a single real condition, namely, an equation for the real parameter b in 
terms of C and D. 

If the map z(() and the constants b, C and D are known, ([51]) gives an explicit expression for F((). Condition (|3T|) 
then provides the following expression for G((): 



G(C) = -ziC^F (C) 



(57) 



This condition is very revealing: it shows that the singularities of G(() dictate the singularities of the conformal 
mapping function z(C). In particular, it shows that if G(C) has a pole singularity at C = then the functional form 
of the mapping function is finite. To see this suppose, for example, that z(£) has the special truncated form 



z(0 



c + 



7 + AO + ^ 



(58) 



fc=i 
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FIG. 4: The geometrical significance of parameters a and c. The fiat interface corresponds to a = 2, c = 0. Increasing 
(decreasing) a above 2 with c = (along the horizontal axis) moves the interface up (down) in a left-right symmetric fashion. 
Non-zero values of c (along the vertical axis) introduce left-right asymmetry. 



where n > 1 is some positive integer. Then 



^ k=i ^ 

Since F(() = 0(C _1 ) as ( ->• 0, it follows from equation (JSTJ) and ([55]) that, near C = 0, G(C) has the form 

G(C) = ^^ + ^ + --- + ^ + Go + G lC + ... (60) 



implying that, near Zd, 



g(z) = - + higher order terms. (61) 



Since our swimmer model means that g(z) must have the form given in (|14|) it is clear that we must pick n = 1. The 
conformal map then has the form 

z(() = -?-+a a + a 1 (, (62) 

C + i 

where a, ao and a\ are complex numbers. 

In addition, there are geometrical constraints on the conformal mapping parameters. Suppose that the singularity 
is at Zd = — i. This implies that z(0) = — i and therefore that 

a = (a — l)i. (63) 

Furthermore, the condition that the interface tends to y — implies 

Re[a x ] = | - 1. (64) 
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Hence, the mapping function can be expressed as 

*(C) 



c + 



T+i(a-l)+(|-l + ic) C, 



(65) 



where a and c are two real parameters. These two parameters have the following geometrical interpretations. When 
a = 2 with c = the interface is flat; increasing a above 2 corresponds to a symmetric upward deformation of the 
interface; decreasing a below 2 leads to a symmetric downward deformation. Changing the value of c away from 
zero introduces left-right asymmetry of the free surface deformation about a vertical axis through the swimmer, as 
illustrated in Figure |H 

The velocity of the swimmer in the co- moving frame is given by the finite part of the expression u + iv at z = — i 
and, for a steady solution, this must vanish. If 

(-is* + d* 



f(z) 



Z + l 



7+/o + /i(z + i) + 



q 



(z + i) 



+ 9a+gi(z + i) + 



the finite part of u + iv at z = — i is given by 

-/o-i7i" + 5T = 
The terms /o, fi and g\ are given (using the residue theorem) as 

no dz 



(66) 



(67) 



fo 



1 



<K, A 



27ri /rO(C) +i) d( 

J_ £ z{l/QF{£) dz 
2?ri 



1 

27ri 



F(Q dz 
(z(C)+i) 2 dC 



9i 



(68) 



d( 



, r (z(C)+i) 2 rfC 

where T is any simple closed curve surrounding f = 0. Note that with determined from equation ([50)l . G(C) 

from (|57|) and /o, /i and from (|68]l . condition (|67]) is a linear equation for constants C and D. It should also be 
pointed out that fo, fi and gi can, in principle, be found in analytical form but the algebra involved is lengthy and 
it is easier to make use of the integral expressions (|68j) to compute these quantities. 

If the values of a, c and U are specified, C and D can be found by simultaneously solving 



-fo- i/i + .9i = 0, 
which is a complex equation, together with the real equation (|55p. i.e., 

J c (-i, Ca) + 2iL> + C + 2Z)i - C = 0, 

and the real equation 



- = Re 
2 



7 cc (-i, Ca) + 6D - 2Ci - 2D 



(69) 



(70) 



(71) 



With constants C and Z? determined in this way, functions F(Q and G(C) (and hence the flow field) are fully 
determined in analytical form. The stresslet, dipole and quadrupole strengths can then be readily computed to be 

2 



D 



q 
d* 



3a - 2 
2 

a - 2 



3a - 2 



(72) 



1C 



LDa 



3 + 5ic - C 



3a - 2 



As a check on the solution scheme, it was verified that when a = 2, c = (so that the interface is flat) we retrieved 
the values 



s* = 2i, q* =0, d* = -4, 



(73) 
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FIG. 5: The two solution branches: the values of c against a are plotted for Ca = 1 and 1000. The upper branches have positive 
values of c, the lower branches have negative values of c. The flat interface parameters are also indicated to show that the 
actual solution branches are disconnected from it. 

which are the values obtained in ^IVI using the method of images. 

We have thus found that there is a two-parameter family of solutions, parametrized by a and c, for a force-free (but 
not necessarily non-rotating) point swimmer - characterized by a singularity description of the form (|14[) - steadily 
translating beneath a deformed free surface with speed U = — 1. 

It remains, however, to enforce the condition Im[/i] = 0. In contrast to the case when the interface is assumed to be 
flat, we now have additional freedoms (in the choice of a and c) to enforce this condition. Physically, this means that 
possible equilibrium solutions exist within a two-parameter family of possible interface shapes. With a two-parameter 
family of possible solutions, and only a single additional requirement, we choose to specify a and see if there are any 
corresponding values of c such that Im[/i] =0. It must also then be checked a posteriori that the map ([650 with 
these values of a and c are one-to-one (univalent) mappings to the fluid region. It has been found that such solutions 
do indeed exist. They are described in the next section. 

VI. CHARACTERIZATION OF THE STEADY SOLUTIONS 

For a given value of a, admissible values of c are found by insisting that 

Im[A] = 0. (74) 

The nature of our formulation is such that the only stage in which we are required to solve a nonlinear equation is 
in satisfying ([740 . Newton's method is used to satisfy this condition by fixing a value of a, taking a guess for c and 
then computing the associated functions F(() and G(£) together with the values of /o,/i and g\. The value of c is 
then updated iteratively until condition (|74[) is satisfied. It is then checked a posteriori that the resulting conformal 
mapping is a one-to-one mapping from the interior of the unit £-disc to the fluid domain. 

Figures [5] shows a graph of c against a when Ca = 1 and Ca = 1000. Focussing on the case Ca = 1 first it is clear that 
there are two solution branches: for certain values of a (but not all) two possible values of c giving physical admissible 
solutions are found. This is not surprising given that a and c are simply mathematically convenient parameters that 
are related to physically meaningful parameters (e.g. the multipole strengths) via highly nonlinear relations. 

The lower solution branch corresponds to c < 0: Figure [5] shows typical profiles throughout the range of existence 
of this branch and reveals that the point of highest curvature is to the right of the swimmer. 

Figure[7]shows the magnitudes of the multipole strengths as a function of a for this branch of solutions. Interestingly, 
as a — > 1.66 which is the lowest a- value for which solutions are found to exist on this branch the free surface profile 
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FIG. 6: Free surface profiles (with Ca = 1) for four different a in the range of existence of solutions for the lower solution 
branch in Figure [S] corresponding to c negative. 
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FIG. 7: Multipole strengths for Ca = 1 in the range of existence of solutions for the lower solution branch with c negative. 

tends to a configuration in which the point of highest curvature tends to a finite value but the multipole strengths 
appear to grow without bound (hence the vertical asymptote shown in Figure [7]). On the other hand, as a tends to 
the upper value for which solutions exist the value of the maximum curvature of the interface appears to grow without 
bound while the values of the multipole strengths appear to asymptote to finite values. 

Figures [5] and |H] show similar graphs for the upper branch for which c is positive. Now the point of highest 
curvature on the interface is to the left of the swimmer position. The qualitative behaviour of the solution branch 
is very similar, however, to the lower branch just investigated. At one end of the range of existence of solutions the 
maximum curvature of the interface tends to a finite value with the multipole strengths growing without bound; at 
the other end of the range, the maximum curvature of the interface grows without bound with the multipole strengths 
levelling off to well-defined values. 

The above results appear to be only weakly dependent on the capillary number Ca. Figure [5] also shows c against 
a in the case Ca = 1000 and once again two distinct solutions branches are found (disconnected from the flat state 
solution). One solution branch has c < 0, a second has c > 0. The free surface profiles for the two solution branches 
are depicted in Figure [10] and are qualitatively similar to the Ca = 1 case. 
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FIG. 8: Free surface profiles (with Ca = 1) for three different a in the range of existence of solutions for the upper solution 
branch in Figure [5] corresponding to c positive. 




a 

FIG. 9: Multipole strengths for Ca = 1 in the range of existence of solutions for the upper solution branch with c positive. 




FIG. 10: Free surface profiles for the two solution branches with Ca = 1000 (left: lower branch; right: upper branch). Shapes 
are qualitatively similar to the Ca = 1 case (Figs. [6] and [8}. 
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It appears therefore that modes of steady translation are possible and that they are associated with a range of 
free surface profiles. What is clear from the solutions is that the steady motion requires significant deformations of 
the surface. In other words, the nonlinear interaction between the swimmer and the dcformability of the interface is 
crucial for steady swimming. 

It is natural to ask which of the family of steady solutions might be realized in practice or be energetically preferable. 
Usually such matters are decided based on some definition of swimmer efficiency. An example is to minimize the energy 
dissipation associated with a displacement of the swimmer by unit distance. Unfortunately, owing to the fact that 
the model used here involves introducing singularities in the flow then the associated dissipation rates are unbounded 
(however, some progress in this direction has recently been made by [37| where use is made of the exact solution 
class found here). Nevertheless, it is intuitively sensible that the particular biological characteristics of any given 
swimmer will mean that there is an upper bound on the strength of the far-field multipoles in any effective singularity 
description of the swimmer. Moreover, it is also sensible to suppose that the rate of energy dissipation will be smaller 
when the multipolc strengths are smaller. Figures [7] and show that the multipole strengths are least when the values 
of a are greatest. As shown in Figures [6] and |8] it is for these higher values of a that the free surface tends towards a 
near-cusped intcrfacial shape (either in front of, or behind, the swimmer). These considerations suggest that efficient 
swimmers in steady motion will be inclined to produce free surface deformations exhibiting a region of high curvature. 



VII. CONCLUSION 

In this paper we have used complex analysis to identify a new mode of low-Reynolds swimming, namely organisms 
or devices that generate flow disturbances which exploit the deformation of a nearby free surface to swim steadily 
parallel to the interface. The resulting surface shapes are illustrated in Figures 161151 and [TU1 The two solution branches 
we obtain, displayed in Figure [5] are seen to be disconnected from the flat state solution. The interaction between 
the swimmer's disturbance to the flow, the free surface, and the subsequent swimmer's dynamics is highly nonlinear, 
and it is likely to be difficult to derive the solutions presented here by a weakly nonlinear perturbation analysis about 
the flat state. It is therefore of mathematical interest that by employing complex variable methods and a singularity 
model we have found these nonlinear interactions to be mathematically tractable. 

The swimmer model we have used is motivated by a simple trcadmilling circular motion of the kind first considered 



by 31] • In principle the mathematical techniques used here arc extendible to more detailed and realistic singularity 
representations. Furthermore it should be possible to extend this study by using the method of matched asymptotic 
expansions to match an "inner" flow generated by a small finite-area swimmer with an "outer" solution in which the 
flow generated by that swimmer interacts with the free surface. The assumptions of such an analysis would be that 
the swimmer size is small compared to its distance from the interface so that, at leading order in an expansion based 
on a parameter expressing this ratio of scales, the swimmer would look like a point singularity. The solutions found 



here would then serve as the outer solutions in such a scheme. 38( has used precisely such a strategy of matched 
asymptotics in his two-dimensional complex variable model of a dcformable bubble in the flow field generated by 
Taylor's four-roller mill. 

The bearing of our two-dimensional analysis on the physics of actual organisms employing a free surface for propul- 



sion is not clear. Other models for free surface locomotion have been proposed [24[ also involving free surface interac- 
tion, but they were restricted to small-amplitude motion. Here, using a fully nonlinear analysis, we have unveiled an 
essential mechanism for the steady locomotion of a swimmer near a free capillary surface. In particular, our analysis 
has demonstrated the essential role played by interface dcformability, and hinted at the significance of regions with 
high surface curvature. 

It is likely that, as in the case of swimming near a no-slip wall, the generic locomotion mechanism of an organism 
near a free surface will have a more complicated spatio-temporal structure. Indeed, free surface deflection associated 
with unsteady undulatory waves propagating along the foot of a water snail have been observed in practice. In future 
work it may be interesting to study the full unsteady dynamics of a point swimmer near a free capillary surface. 
The geometrical complexities even in that simplified case would no doubt require a fully numerical investigation. An 
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analysis of the stability of the steadily translating equilibria just identified is also of interest but is left for a future 
study. 
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